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We derive the gravitonic Casimir effect with non-idealised boundary conditions. This allows the 
quantification of the gravitonic contribution to the Casimir effect from real bodies. We quantify 
the meagreness of the gravitonic Casimir effect in ordinary matter. We also quantify the enhanced 
effect produced by the speculated Heisenberg-Couloumb (H-C) effect in superconductors, thereby 
providing a test for the validity of the H-C theory, and consequently the existence of gravitons. 
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One of the most remarkable consequences of the non¬ 
zero vacuum energy predicted by quantum field theory, 
is the Casimir effect. In its most basic form, the Casimir 
effect is the attraction between two perfectly reflecting 
surfaces as a result of the restriction of allowed modes 
in the vacuum between them (Fig. [[]). Real bodies how¬ 
ever are not perfectly reflecting, and the generalisation 
of these ideal boundary conditions to more realistic ones 
have been derived for the electromagnetic (EM) field, re¬ 
sulting in the Lifshitz formula at zero temperature [lj. 
The EM field of course, is not the only field that pro¬ 
duces the Casimir effect; in theory all fields of the quan¬ 
tum vacuum contribute to the Casimir effect. In fact 
the contribution to the Casimir effect from any massless 
field which is opaque to the plates should be significant 
(mass quickly weakens the Casimir effect !)■ Therefore 
one may imagine plates which are opaque to the gravi¬ 
tational field, so that the Casimir effect would then be 
a manifestation of the quantisation of the gravitational 
field, or gravitons. The difficulty is in finding such a 
medium, as ordinarily materials are transparent to the 
gravitational field Q ■ 


Recently however, there have been suggestions that the 
properties of quantum fluids (superconductors, superflu¬ 
ids, quantum Hall fluids, Bose-Einstein condensates) may 
enhance the interaction with gravitational waves (GW). 
The novel effects of the interaction of a gravitational field 
with a quantum fluid was first investigated by DeWitt [|] 
and Papini E, who calculated that a Lense-Thirring field 
should induce a current in the superconductor. Following 
this, further analyses were made into the interaction of 
GW with superconductors EQ, proposing superfluids 
as a medium for gravitational antennae |8|], supercon¬ 
ducting circuits as GW detectors @, transducers Him 
and mirrors 12]. These idea have not been met without 
controversy [l3l . El- Although a few experiments have 
attempted^ test the proposed enhanced GW interac¬ 
tion 15. .EH, none have produced clear and unambiguous 
outcomes. This is perhaps because of the small magni¬ 
tude of some of the theorised effects coupled with the 
practical challenges in producing the environment capa¬ 
ble of their detection jl6(. In light of the controversial 
status of enhanced gravitational interaction, an under- 



FIG. 1. The Casimir effect with a two plate setup. The change 
in the refractive index of the plates causes the gravitational 
wave to refract, k represents the wave vector of the incident, 
transmitted, and reflected gravitational waves, and 7 is the 
corresponding angle with respect to the surface normal. 


standing of the gravitonic Casimir effect for realistic bod¬ 
ies is needed to accurately quantify the contribution in 
ordinary materials and in any theorised system with en¬ 
hanced gravitational interactions. 


The Casimir effect has also been investigated in weak 
gravitational fields to see the effect the slightly curved 
spacetime background would have on the Casimir en¬ 
ergy 17]. In these works, the gravitational field is as¬ 
sumed to be unaffected by the matter that form the 
boundary conditions. This is different from what we 
are considering here, where we look at how the gravi¬ 
tational field of the vacuum interacts with the matter 
of the boundary conditions to give rise to the gravitonic 
Casimir effect, in flat spacetime. This gravitonic Casimir 
effect has been considered in a cosmological context; how¬ 
ever in these works, the boundary conditions are idealised 
and not suitable for realistic terrestrial systems []j| . The 
gravitonic Casimir potential also has been calculated for 
a massive test point particle interacting with a fluctu- 
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ating mass distribution [19j. In this letter we derive the 
gravitonic Casimir effect for real bodies. In effect we give 
the Lifshitz formula for the gravitonic contribution to the 
Casimir effect at zero temperature. 

The contribution of small perturbations to the flat 
spacetime metric, 77 ^, is well described by the linearised 
Einstein field equations. We begin with a Maxwell¬ 
like formulation of the linearised Einstein field equations 
known as gravitoelectromagnetism (GEM) 0 E], 


= Kp {E) , 

( 1 ) 

= «P (M) , 

( 2 ) 

-— - «J (M ) 
dt ’ 

(3) 

<9E -r(F\ 
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where k = 8nG/c 4 . In this formulation, components of 
the Weyl tensor (C a ptiv) play roles analogous to the elec¬ 
tric and magnetic fields in electromagnetism: Eij = CVb.rw 
and Bij = * 60 * 07 . ( * denotes Hodge dualisation 122l| V 
The rhs of the GEM equations contain components of 
the matter current, J pvp = -T p ^^ + where 

T pv (T = Tjf) are the stresses due to the perturba¬ 
tion. /?* = A 00 , Pi = * -/,oo, Jij — JiOj , and 

= *Jioj (we use the convention that Greek letters 
go from 0 to 3 and Roman letters go from 1 to 3, unless 
otherwise stated). 

We are interested in the macroscopic properties of the 
system, so E, B, T will represent Russakoff spatial aver¬ 
aged values [23j, as is used in the macroscopic Maxwell’s 
equations. We assume that T Z j = \{u)Eij, where \(lo) is 
the frequency dependent gravitational susceptibility [24|; 
we will later show a specific model of matter where 
the induced stresses are of this form. We consider the 
gravitoelectromagnetic fields to be of plane wave form, 
E i:j (r,t) = £i j e i ^' r ~ ult \ B i3 (v,t) = B i3 e i ( k ‘ r -“ t ). From 
the GEM equations, E and B are transverse waves with 
two independent polarisations, “+’ and ‘x’, respectively 
defined by the only non-vanishing components £\ \ = 
—£■22 = B 12 = B 21 = o. and —Bn = B 22 = £12 = £21 = 
ft, in the proper frame of the plane wave. Here we define 
the proper frame of the plane wave as that in which k is 
along the positive 2 -axis. 

The vacuum energy of any quantum field between par¬ 
allel plates (separated by distance a) is a summation of 
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the energy of all allowed modes of the field. The opaque 
boundary of the parallel plates in the sy-plane means 
that the k z components of the field are discrete, whereas 
the k|| = (k x ,k y ) components remain continuous. In 
terms of the graviton eigenfrequencies, the vacuum en¬ 
ergy of gravitons between the plates at zero temperature 
is given by [25J, 




k \\ dk \\ + UJ n) cr , ( 5 ) 


where a is the surface area of the plate. The allowed 
modes (w+, w*) between the plates are found by consid¬ 
ering the boundary conditions, as follows. 

A naive application of Stokes’ theorem to the E and 
B fields at the plate interface will produce an overde¬ 
termined problem. The reason for this is that the extra 
components of these second-order tensors would intro¬ 
duce extra constraints, as compared to vector fields, such 
as the EM field. Thus E and B cannot be considered 
completely smooth across the interface. Instead only the 
traceless part of the tangential components of the tensor 
fields are considered smooth across the interface - this is 


known as the smoothness principle 24). This gives rise 
to the boundary conditions, 


A[(l + KX /2)E il ]=0 
AB tt = 0 , 


( 6 ) 

(7) 


where E ee = £?t e i(k r-ujt) - g (-j le traceless part of E 
after it has been projected onto the interface, and B// e = 
Bab'e dk is the traceless part of B after it has been 
projected onto the interface (subscripts a,b = 1,2). In 
the interface frame as shown in Fig. (HI), £-j r 1 T = —£ 2 2 = 
a(l — S 2 / 2), £ E 2 = £' 21 ' = PC, and Bf E = —Bj/f = 
—/3(1 - S 2 / 2 ), B{ 2 t =H 2 T = aC, where S = sin 7 = 
kJk and C = cos 7 = k z /k. A Q = Q 1 — Q 2 refers to the 
change in quantity Q at the interface between medium 1 
and medium 2. 

We now adapt van Kampen et aUs (26j contour in¬ 
tegral method, except with different boundary condi¬ 
tions, to get the gravitonic Casimir energy. Applying the 
boundary conditions of Eq. © and 0 at the two plate 
interfaces, we get the following system of linear homo¬ 
geneous equations of variables a, a',a", 0 !" for the 
polarisation (primes indicate the region of operation, as 
represented in Fig. |T|), 


a'(l + «x/2)(l - S' 2 /2)e~ qa / 2 = a(l - S 2 /2)e~ qa / 2 + a" {l - S 2 /2)e qa/2 , ( 8 ) 

a'"(l + k X /2)(1 - S' 2 /2)e~ q ' a/2 = a(l - S 2 /2)e qa/2 + a "{\ - S 2 /2)e~ qa/2 , (9) 

a'C'e~ q ' a/2 = aCe~ qa/2 - a"Ce qa/2 , (10) 

- a'"C'e- q ' a I 2 = aCe qa / 2 - a"Ce~ qa/2 , (11) 
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where q 2 = —k 2 = fc 2 — w 2 /c 2 , and q' 2 = —k' 2 = /cjj — 
(1 + k X )uj 2 /c 2 . A non-trivial solution exists when the 
determinant of the corresponding matrix is zero, 


The boundary conditions Eq. © and © at an in¬ 
terface yield Eq. © and (in. Simultaneously solving 
Eq. © and (flOl) one gets 


A+ = e~ aq '{[C'(S 2 - 2) - (1 + k X /2)C(S' 2 - 2)} 2 e~ aq 

- [C'(S 2 - 2) + (1 + k X /2)C{S' 2 - 2)] 2 e° 9 } = 0 . (12) 

The boundary conditions also yield a similar set of linear 
homogeneous equations of variables /3,/3',/3from 
which we get for the 1 x ’-polarisation, 

A x = e" a9 '{[(l + k X /2)C'(S 2 - 2) - C{S' 2 - 2)} 2 e~ qa 

- [(1 + k X /2)C\S 2 - 2) + C(S' 2 - 2)]V a } = 0 . (13) 

Solutions of Eq. m and m give the allowed modes 
between the plates, i.e. the eigenfrequencies (w+, w x ). 
There will be an infinite number of these eigenfrequen¬ 
cies. We can sum them in Eq. ©, using the argument 
principle of complex analysis [271 ]. 


E 




n 


1 

27 ii 


—zoo 

J uid In A + 


wdln A 


_ 200 


C+ 


(14) 


where the closed contour of integration has been taken 
over a semicircle C + of infinite radius in the right half 
of the complex plane of to with its centre at the origin, 
and the imaginary axis [too, —zoo], in a counterclockwise 
manner. 

In the high frequency limit, the system does not 
have time to respond to the rapid oscillations of the 
field, and therefore will effectively act with time av¬ 
eraged behaviour. Analogous to the EM case, we 
thus make the natural assumption lim^^oo X (cd) = 
liniu^oo d X (ui)/duj = 0 and therefore 7 = 7 ' in this limit 
i.e. at very high frequencies the plates are effectively 
transparent to the GW. With this assumption, the sec¬ 
ond integral in Eq. m is independent of a. 

The summation over the infinite number of allowed 
modes will of course give an infinite value for Eq. To get 
the finite Casimir energy per unit area we take away the 
energy at infinite separation. Employing the argument 
principle to Eq. © and then integrating by parts, one 
arrives at the gravitonic Casimir energy (£ = — ico), 


E{a) 


Eq Eq 

-lim — 

(7 a—>00 (J 

h f°° 

k\\dk\\ 


47r 2 


r o 


f‘00 

/ [ln(l - r 2 e~ 2qa ) 

Jo 

+ ln(l -r 2 x e- 2qa M, 


(15) 


where, 

_ C(S 2 - 2) - (1 + k\/2)C{S' 2 - 2) 
r+ = C'(S 2 - 2) + (1 + K X /2)C(S r2 - 2) ’ 
_ (1 + k X /2)C'(S 2 - 2) - C\S' 2 - 2) 
rx ~ (1 + kx/2)C’(S 2 - 2 ) + C(S' 2 - 2 ) ' 


(16) 

( 17 ) 


_ -CJS 2 - 2) + (1 + K X /2)C(S' 2 - 2) qa 
a C'{S 2 - 2) + (1 + n X /2)C(S' 2 - 2) ' 1 ’ 

Comparison of Eq. (fl 8 l) with Eq. HU shows that, |r + | = 
|a"/a|, i.e. r + is the magnitude of the gravitational 
reflection coefficient of the ‘+’-polarisation. Similarly 
for the ‘ x’-polarisation, |r x | = \fd"/fd\ is the magni¬ 
tude of the gravitational reflection coefficient of the ‘ x ’- 
polarisation. This tells us that Eq. m and m are 
the Fresnel reflection coefficients for the two polarisations 
of the gravitational wave. Therefore Eq. m has the 
same form as the Lifshitz formula for the EM Casimir 
energy at zero temperature, except the EM reflection 
coefficients have been replaced with their gravitational 
equivalent. It is important to point out that this was 
not a priori obvious, as our starting point was the lin¬ 
earised Einstein equations, which is fundamentally dif¬ 
ferent from Maxwell’s equation of electromagnetism, al¬ 
though GEM provides useful analogies between the two 
theories. A marked difference is that the EM fields are 
vectors whereas E and B are tensor fields. 

We need now to find an estimate for ,\(w). A number of 
models of the interaction of GW with bulk matter exists. 


For example, Ref. [2lj,[28] considered a medium consisting 
of molecules modelled as individual harmonic oscillators 
to calculate the quadrupole moment induced by an inci¬ 
dent GW; Ref. [29( studied the interaction and dispersion 
of GW in a hot gas. Here we use Peters’ [30] model who 
considered the scattering of GW by the gravitational field 
of individual free particles of a thin sheet of matter. Pe¬ 
ters derive a gravitational refractive index n which was 
much larger than that generated by just considering the 
induced quadrupole moments, suggesting that his model 
encapsulates the dominant GW interaction with matter. 
Peters gives the gravitational refractive index as, 


n= 1 + 


27 xGp 


(19) 


where p is the density of the medium. 

Under the Lorentz gauge (/i M “ jCt = 0), d a d a h^ v = 
—2 kT^ (note is traceless). Using this, Eq. (fl9l) . 
and the relation Eij = —uj 2 hij/2 , result in gravitational 
stresses of the previously assumed form, = X (co)Eij 
where, 


xM = 


1 — n{uj) 2 

KC 2 


( 20 ) 


with the high frequency limits: lim^^oo X (co) = 
lim^oo d X (uj)/du! = 0 . 

Except for the lowest frequencies, for ordinary mate¬ 
rial parameters, X <C c 4 /G, and therefore the reflection 
coefficients and Casimir pressure [P(a) = — dE(a)/da\ 
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are negligible. Specifically, if we take a typical ()\p] = 
10 4 kg/m 3 and 0[l\ = lA, then P(1CT 6 ) « 10" 21 nPa. 
This value is of course beyond detection. For the Casimir 
pressure to be measurable, the density of the material 
would have to be at the very least 0\p] = 0{c 2 /G\ ~ 
10 2 'kg/m 3 . This material density is clearly not achiev¬ 
able, at least terrestrially. 

Up until now we have only considered materials with 
classical properties. Recently, Minter et al. |12] pro¬ 
pose that the quantum mechanical properties of super¬ 
conducting films can give rise to specular reflection of 
GW. Their claim is that in an ordinary metal plate, the 
ions and normal electrons locally co-move together along 
the same geodesics in the presence of a GW. However, 
when the plate becomes superconducting, the quantum- 
mechanical non-localisability of the negatively charged 
Cooper pair undergoes non-geodesic motion, whereas 
the positive charged ions of the lattice remains on the 
geodesic path. Specifically, the authors couple the grav¬ 
itational field to the superconductor through DeWitt’s 
minimal coupling scheme jSJ. To first-order, the ground 
state of the delocalised Cooper pairs do not change in the 
presence of a GW whose frequency is less than the BCS 
gap frequency. In comparison, the shift in the momentum 
of the localised ions is proportional to the gravitational 
vector potential. Because the Cooper pairs and ions are 
oppositely charged, a strong Coulomb force will resist this 
separation of charge caused by the GW, resulting in its 
reflection. The authors dub this the Heisenberg-Coulomb 
effect. A similar effect had been proposed in Ref. jfi}. 

The origins of the arguments employed by Minter et 
al. are heuristical in nature, some of which we believe 
require a much more formal approach to be convinc¬ 
ing. This is echoed in a review article 0 on theories 
of enhanced gravitational interaction with quantum flu¬ 
ids, which predates Ref. [12]; in particular the authors 
urge that a formal derivation of the quantum mechanical 
coupling between an electron and both the electromag¬ 
netic and gravitational field is needed. Nevertheless, the 
work by Minter et al. do yield results which can be used 
to falsify their theory. The H-C effect should enhance the 
Casimir pressure between superconducting plates. Here 
we quantify the size of this effect. 

Minter et al. give the reflection coefficient of a super¬ 
conducting film from an incident GW as, 




( 21 ) 


where <5 is the EM skin depth of the superconducting 
film and d the film thickness. In the thin superconduct¬ 
ing film, the current flows in the x — y plane; therefore 
only the normal incident component of the GW will drive 
the current. At normal incidence the magnitude of the 
reflection coefficient of the two polarisation modes are 
the same, as they only differ by a rotation (this is also 


true for EM waves). We can use Eq. m in Eq. USD 
to calculate the gravitonic Casimir pressure between two 
superconducting films. 

We calculate the gravitonic contribution to Casimir 
pressure for superconducting lead (Pb) of thickness d = 2 
nm at zero temperature. The EM skin depth of Pb is <5 = 
37 nm. We compare this with the photonic contribution 
to the Casimir pressure of superconducting lead. The 
EM reflection coefficient is 0 


te 



( 22 ) 


where A = 83 nm is the coherence length. The photonic 
contribution to the Casimir pressure is calculated by us¬ 
ing Eq. (B^l) in the EM Lifshitz formula [l|, which has 
the same form as Eq. m- 

Eq. (l22l) and (I2T1) are most valid when the driving fre¬ 
quency is less than the BCS gap frequency, as for higher 
frequencies the dissipative component of the complex 
conductivities would need to be taken into account. How¬ 
ever, as higher frequency modes contribute exponentially 
less to the Casimir energy, one may still use the simple re¬ 
flection coefficients derived by Minter et al. to obtain the 
first-order estimate of the Casimir pressure for supercon¬ 
ductors at zero temperature. Fig. [2] compares the gravi¬ 
tonic to photonic contribution of the Casimir pressure 
as a function of separation of plates of superconducting 
Pb. It shows that gravitons can have a significant con¬ 
tribution to the Casimir pressure, via the H-C effect. In 
fact, for the superconducting Pb film considered here, the 
gravitons will dominate the Casimir pressure by an order 
of magnitude over photons. 

The magnitude of the Casimir pressure and plate sepa¬ 
ration distance that we are talking about here, is compa¬ 
rable in size to what has already been achieved in current 


experiments [31l 1321 ]. Few experiments however, have 
been conducted at low temperatures j33|, as room tem¬ 
perature setups are a more experimentally accessible en¬ 
vironment. Of these low temperature investigations, only 
the ALADIN project has experimented with supercon¬ 
ducting aluminium (Al) film, separated by a thin oxide 
layer from a thick gold plate, to observe how the Casimir 
energy influences the superconducting phase transition 
(preliminary experimental results are reported in 34jf. 
One could imagine modifying such an experiment to test 
for the gravitonic contribution to the Casimir effect as 
described in this letter. 

We have derived here a Lifshitz-type formula for the 
gravitonic Casimir effect for real bodies. Besides com¬ 
pleting a theoretical gap in our understanding of the 
Casimir effect, this formula is important in light of re¬ 
cent models of enhanced gravitational interaction, as it 
allows us to quantify the gravitonic contribution to the 
Casimir effect predicted by these theories. If measure¬ 
ments of the Casimir pressure of the setup described in 
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a [pm] 

FIG. 2. Gravitonic (dotted line) and photonic (solid line) 
contributions to the Casimir pressure of parallel plates of 
superconducting Pb at zero temperature, as a function of 
plate separation a. Due to the Heisenberg-Coulomb effect, 
the gravitonic contribution exceeds the photonic. 

this letter match only to a photonic contribution [Fig. © 
solid line], then one should conclude that the H-C effect 
is invalid, if we are to hold on to the idea of the graviton. 
However, if experiments show the Casimir pressure to be 
an order of magnitude larger than that predicted from 
the photonic contribution alone, this would be the first 
experimental evidence for the validity of the H-C theory 
and the existence of gravitons. This would open a new 
field in the way of graviton detection. 
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